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Based on the SU(2) lattice gauge theory formulation of the t-J model, we discuss possible signature 
of the unit cell doubling associated with the staggered flux (SF) state in the lightly doped spin liquid. 
Although the SF state appears only dynamically in a uniform d-wave superconducting (SC) state, 
CO ' a topological defect [SU(2) vortex] freezes the SF state inside the vortex core. Consequently, the 

unit cell doubling shows up in the hopping (xij) ^-nd pairing (Aij) order parameters of physical 
electrons. We find that whereas the center in the vortex core is a SF state, as one moves away from 
the core center, a correlated staggered modulation of Xij a-nd Aij becomes predominant. We predict 
, that over the region outside the core and inside the internal gauge field penetration depth around a 

^ ' vortex center, the local density-of-states (LDOS) exhibits staggered peak-dip (SPD) structure inside 

the V-shaped profile when measured on the bonds. The SPD structure has its direct origin in the 
. unit cell doubling associated with the SF core and the robust topological texture, which has little 

to do with the symmetry of the d-wave order parameter. Therefore the structure may survive the 
tunneling matrix element effects and easily be detected by STM experiment. 

o 

o . 

; I. INTRODUCTION 

(N 

t> ■ High Tc superconductors are doped Mott insulators. Soon after the discovery, Anderson proposed that the strong 

. correlation physics of the doped Mott insulator is well captured by the t-J model. Taking account of competition 

, between the hole kinetic energy xt and the spin exchange energy J, he proposed that the spin liquid states formed 
out of the resonating valence bond (RVB) singlets are a good starting point to study this model. tl A standard way 
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of enforcing the constraint of no double occupancy in the t-J model is the slave boson formalism where a physical 
' electron operator Cia with spin a at the site i is splintered into auxiliary spin- 1/2 fermion fi„ and charge- 1 boson hi. 
' Cia — fiabl- One hue to describe the spin liquid state is to start from mean field (MF) decoupling, i~i Aij = (eaafiafia) 
^ ■ and Xij = ifjafic') which characterize the spin liquid state formed out of the RVB singlets. The phases of Xij ^-^d 
• Aij transform as the lattice gauge fields under local U(l) transformation which naturally leads us to a U(l) gauge 
5 ; theory.!! At zero doping, the t-J model reduces to an antiferromagnetic Heisenberg model which has an exact local 
^ SU(2) gauge symmetry.! Then, the translationally invariant solution can be described as a 7r-flux state^l or a d-wave 
^ \ pairing statei with \xij\ = These apparently different mean-field ansatz describe exactly the same MF state, 

[ since they are just SU(2) gauge equivalent. 
$H In the U(l) slave-boson formulation, however, the SU(2) symmetry is broken upon hole doping due to appearance 

■ - - ' of the boson hopping term. Consequently, the d-wave superconducting (SC) state and the flux state are no longer 
equivalent. For small doping and small J/t, the 7r-flux phase at zero doping is disfavored against the staggered 
flux (SF) phase with \xij\ > |Aij|.i~0 The SF state, however, breaks physical symmetries associated with the time 
reversal and the spatial translation, which causes the unit cell doubling and staggered orbital currents of the physical 
holes. Eventually the SC phase is picked out as the MF solution out of a inflnite number of degenerate states upon 
doping. However, it is still quite natural to expect the SF state is nearly degenerate with the SC state in the 
lightly doped spin liquid states. That is to say, as far as we conflne ourselves to the spin liquid state, the SU(2) gauge 
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structure at zero doping may still be useful to describe the low energy states in the underdoped regime, which are 
missing in the U(l) formulation. To substantiate this idea, Wen and LeeEllli introduced an SU(2) boson doublet 
{hi)-^ = {bii,bi2) and constructed the effective model which recovers local SU(2) symmetry even upon doping. From 
this viewpoint, the SF state plays a crucial role to describe the low-energy spectrum of the lightly doped spin liquid 
state. The question which we must consider next is how to detect a signature of the SF state contained in the low- 
energy excitation spectrum. The first step in this direction was addressed by Ivanov, Lee, and Wen0 who found a 
signature of the staggered current-current correlation by using a Gutzwiller-projected c?-wave pairing wavefunction. 
This is naturally interpreted as a consequence of the quantum fluctuations around the SC state toward the SF state. 
LeungEl further sought for a signature of the SF state and found the current-current correlation in the d-wave SC 
state by using exact diagonalization of the t-J model for a system with two holes on a 32-site lattice. 

In the experimental side, structure of the low-energy excitations in the underlying "normal" metallic phase is con- 
cealed by a phase transition to bulk superconductivity. One promising way to escape from this situation is to introduce 
the topological defect into the superconducting phase, i.e., the vortex. Inside the vortex core, low-energy properties 
of the normal metallic phase show up against the surrounding superconducting phase. Remarkable progress in low- 
temperature STM technique with atomic resolutionEl has given us good opportunities to look into the electronic 
states around the superconducting vortex. 0^0 Recent STM experiments on Bi2Sr2CaCu208 (BSCCO)0@ revealed 
the striking fact that the normal core electronic state exhibits the "pseudogap" structure characteristic of the normal 
state pseudogap above Tc- A description of a vortex core based on conventional BCS theory requires that the super- 
conducting order parameter vanishes inside the core, which is usually accompanied by the vanishing of the energy 
gap. The experimental finding thus strongly suggests that the electronic structure of the vortex core is qualitatively 
different from that given by conventional picture. 

The theoretical description of the normal core in the light of the strong correlation physics, however, remains 
unresolved.&@ In the SU(2) picture, since the SF state is nearly degenerate with the SC state, it is naturally 
expected that by frustrating the SC state the SF state will be revealed inside the core. Based on this idea, Lee and 
Wei]0 proposed a model of the vortex with a SF core, characterized by a pseudo gap and staggered orbital current. 
Quite recently, Han, Wang, and Lee found evidence of the SF order near the vortex core by using Gutzwiller projected 
U(l) slave-boson mean-field wave function. These numerical results so faroEjCJ'tJ strongly suggest that the SF 
state is a key ingredient in the t-J model. 

The vortex with the SF core [SU(2) vortex] offers us an opportunity to experimentally detect the SF state at low 
temperatures below Tc, whereas it may be difficult to probe the staggered current pattern in the zero- field uniform 
SC state because of spatial and temporal fluctuations. Possible experimental tests of the SF core were proposed as 
summarized below. (1) Cyclotron resonance or Shubnikov-de Haas experiments in a high quality underdoped sample 
aX H > Hc2 can detect the small Fermi pockets around (±7r/2, ±7r/2) points with not uniformly spaced Landau levels. 
(2) /i-SR or neutron scattering experiments can directly detect the staggered currents which produce a small staggered 
magnetic field of order 10 gauss.S Intensity of the signal may increase upon increasing iJ, since the increasing H 
excites more vortices with the core size being independent of H . (3) NMR experiments can detect side bands in the Y 
NMR line in YBCO samples with a splitting independent of H but with weight proportional to H . For this purpose, 
Y2Ba4Cu70i5 may be ideal, because there are asymmetric bi-layers where the Y ion sits in between, and it may be 
possible to have one plane of the bi-layer optimally doped while the other plane (next to the double chain) remains 
underdoped, i.e., the staggered magnetic field at the Y site does not cancel. 

Now we are naturally lead to the following question: is it possible to detect a signature of the unit cell doubling 
associated with the SF core through the state-of-the-art STM technique? It turned out that there is no effect inside 
the SF core, because what is staggering in the SF state is the currents, which does not show up in the charge density. 
This situation motivated us to look at the region outside the core. We addressed this problem in our previous paperEa 
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and found that whereas the center in the vortex core is a SF state, as one moves away from the core center, a 
correlated staggered modulation of the hopping amplitude Xij ^-nd pairing amplitude Ay of the physical electrons 
becomes predominant. We predicted that in this region, the LDOS exhibits staggered modulation when measured on 
the bonds, which may be directly detected by STM experiments. 

In this paper, we give a full account of the results summarized in Ref. 1^1] and examine the LDOS around the SU(2) 
vortex in detail. The outline is as follows. In Sec. II, we will give an overview of the SU(2) lattice gauge theory 
formulation of the t-J model (Sec. II A) and then discuss the topological texture of the SU(2) boson condensate 
based on the 0(4) cr-model (Sec. II B). We are mainly concerned with the LDOS outside the core through which we 
detect the unit cell doubling stabilized by the robust topological texture. For this purpose, a close study of the vortex 
core state is not necessary. To take account of the phase winding, we will apply a simple London model for a single 
vortex to the SU(2) vortex model (Sec. II C). In Sec. Ill, we discuss the hopping and pairing order parameters of the 
physical electron around the vortex. For this purpose, we perform an appropriate local SU(2) gauge transformation 
(Sec. Ill A). Then, we argue in detail that as one moves away from the core center, a correlated staggered modulation 
of Xij and Aij becomes predominant (Sec. Ill B). In Sec. IV, we evaluate the LDOS outside the core. Formulation of 
the LDOS at an arbitrary point on lattice is given in Sec. IV A. It is demonstrated that the LDOS exhibits conspicuous 
staggered pattern only when measured on the bonds. To obtain the LDOS, we compute the lattice propagator by 
using two complementary approaches, which are presented in Sec. IV B and Sec. IV C. Finally, concluding remarks 
are given in Sec. V. 



II. SU(2) VORTEX WITH THE STAGGERED FLUX CORE 



In this section, we recapitulate the SU(2) lattice gauge theory formulation of the t-J model and then discuss the 
SU(2) vortex model in some detail. 



A. SU(2) lattice gauge theory formulation of the t-J model 



The t-J model Hamiltonian is given by 

^ = E + H-C.) + J E (^^- - (2-1) 

<i,j>,cr <i,j> 

where and Ci^ are the projected electron operators with the constraint rii < 1. In the SU(2) slave-boson 

approach,EBEl a physical electron is represented as an SU(2) singlet formed out of the "isospin" doublets of the 
fermion (ipi^) and boson (hi): 

= -^hji,,, - 71 (^^^1^- + ^^'^bljl), (2.2) 



^-=(/^^.tV (2-3) 



with 

The physical hole density (&J]^&ii + &j2^i2) = a; is enforced by the chemical potential ^. We need to introduce the 
temporal component of the gauge field aoi to ensure the projection of the Hilbert space onto the SU(2) singlet 
subspace, {^i/jj^T^ia + /ijT/ii)|phys>= 0, which is identical to that of the original t-J model. The conventional U(l) 
slave-boson bi is now regarded as the SU(2) boson doublet having only its isospin "up" component: (h'f^)'^ — (6^,0). 
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The spin liquid state is characterized by the order parameters Ay = {eaafiafjs) and Xij = {fjafjr^) which constitute 



a 2 X 2 matrix 



U^, = ( ] ■ (2.4) 



We should 



By this decoupling, the spin-exchange term is replaced with Si ■ Sj ^ ^ J2a ^ia^ij''t'j<^ + Tf'^'^ 
stress here that in the presence of the &2-boson Xij and Ay cannot be interpreted as the hopping and pairing order 
parameters of a physical electron [see Eq. (2.2)]. The "phase" of Uij is now interpreted as the SU(2) lattice gauge 



fieldsffl 

tJij = Uij exp[-iajj • -r], (2.5) 

where t = {t^ , r^, t^) are Pauli matrices and = (aj^-, of^, af^) is the gauge field on every link. Now the t-J model is 
described by the fermion-boson system interacting with the SU(2) lattice gauge field, described by the Lagrangian 
■.L^ = Ll + LE + i E<.,,> Tr[t7;^.i7„] with 



L^ = ^h\ [5y (a, - /i) + tU,j\ + X hliao^ ■ Th„ (2.7) 

where J = 3 J/8 and t — t/2. The mean- field solution is obtained by integrating out the fermions and minimizing the 
mean-field energy E{{Uij, hi}), which leads to Uij on the links and the boson hi on the sites. 

The SU(2) gauge invariance is realized through the relation E{{Uij,hi}) — E{{WiiJijWj ,Wihi}), for any Wi G 
SU(2). Thanks to the SU(2) symmetry, we can choose a convenient gauge fixing to describe the MF state in an SU(2) 
invariant way. Convenient gauge choices in the underdoped regime are the "d-wave gauge" or the "staggered- flux (SF) 
gauge" specified by 

Uf^=~Xor' + i-iy-+'-A„T\ (2.8) 
C/f/ = -At3 exp[i(-l)*-+-'''$oT^], (2.9) 



respectively, where A — \/ Xo ~^ and $o — tan~^(Ao/xo)- Eq. (2.8) describes fermions with d-wave pairing order 
parameters, while Eq. ( ^.g| ) describes fermions hopping with flux ±4$o on alternating plaquettes.i At zero doping 
(x = 0) there is no boson and these apparently different mean-field ansaz describe exactly the same MF state, since 
Uf^ and t/f/ are just SU(2) gauge equivalent, i.e., Uf^ = w^Uffw] [E{{Uf^}) = E{{w^Uffwl}) ^ E{{Uff})], where 
the transformation is explicitly given by 

^ exp[i{-iy'^+'y^T^]. (2.10) 

Upon doping, however, the Ufj with the U(l) boson condensate (/loi)"^ = {bi,0) — (v^, 0) characterizes the physical 
d-wave SC state, while the ^7,^^ with the U(l) boson condensate (/loi)"^ = (v^j 0) characterizes the physical SF state. 
These states are no longer physically equivalent because of the presence of the boson condensate [E{{Ufj,hQi}) ^ 
E{{lJff ,hQiY)\ and the SC phase is picked out as the MF solution.i0 Accordingly, the "flux"<i>o = tan"i(Ao/xo) 
decreases from $o = 7''/4 (7r-flux phase) upon doping.i0 

The advantage of the SF gauge is that it is apparent that the SU(2) symmetry has been broken down to the residual 
U(l), which we denote U(l)rGs, since Uff contains only r^.Eil The lattice gauge fields and af^- become massive 
by the Anderson-Higgs mechanism and can be ignored, while af^- remains massless and is the important low energy 
degrees of freedom which should be included, i.e. we consider 
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Uff = -Ar^exp [i{^iy-+^y ^qt"^] exp [-mf^-r^] . (2.11) 

In this gauge, we can discuss a vortex structure under the external magnetic field in a way quite similar to the 
conventional BCS vortex where the gauge structure is characterized by only the electromagnetic (EM) U(l)om- The 
difference is that, in our problem, the gauge structure is characterized by U(l)cm®U(l)ics- 

B. 0(4) (T-model description of the local boson condensate 

In the presence of a magnetic field, the mean-field solution contains vortices. The SU(2) vortex mode]0 was 
discussed based on the 0(4) (T-model description for a slowly varying boson condensate. El The basic idea is that at 
low temperatures the bosons are nearly condensed to the bottom of the band and are slowly varying in space and time. 



The ansatzs (2^) and ( ^.9| ) gives the one-boson dispersion = —iA{cos'^ + cos^ ky + 2 cos 2$o cos fc^; cos fcj,)^/^. 
The bi and &2 bosons are then nearly condensed to the band bottom (0,0) and (7r,7r), respectively.^ On the other 
hand, the fermions are fluctuating over the lattice scale and can be integrated out, after choosing an a^i field which 
minimize the action locally. This view is in the spirit of the Born-Oppenheimer approximationJll In the SF gauge 



given by Eq. (2.9), the local boson condensate (LBC) can be written as 



where Zn and z;2 [CP^ fields] are slowly varying in space and time and parameterized by 

2,1 =e^'""cos|, z,2=e^'^-sin|, (2.13) 

with the internal phases being given by 

Vii = tti - ip2i = o-i + 0i/2. (2.14) 

We shall give some remarks on the expression Eq. ( p.l2[ ) in appendix A. 

The overall phase angle a is associated with the U(l)cm- The internal SU(2) gauge symmetry is broken down 
to U(l)res and the angles ^ and d are interpreted as polar angles of the manifold of the LBC: SU(2)/U(l)rcs — S^. 
Topological stability of vortex formation is indicated by the non-trivial topology, 7r2[SU(2)/U(l)rcs] ~ tti [U(l)rcs] — Z. 
The internal degrees of freedom of the LBC is visualized by the vector 

li = zjTZi — (sin 0i cos sin sin COS ^i), (2-15) 

which has the meaning of the quantization axis for the z fields, {zi)'^ = {zn, 2^2). In the SF gauge, the uniform d-wave 
SC state and the uniform SF state are described by 9i — tt/2 and 6i — 0, tt, respectively. The angle 0i is associated 
with the residual gauge symmetry U(l)res which is further broken down to {0} upon the bose condensation which 
triggers the superconducting phase transition. 

The low energy dynamics of the LBC is described by an anisotropic 0(4) cr-model coupled to the gauge fields.El 
Since we are only concerned with static configurations, we shall ignore the time dependent terms from now on. The 
free energy associated with this model is written in a form Fcff = Fk + F± + Fa + Fa, explained below. In the SU(2) 



formulation, only the boson can carry charge. Under the magnetic field, the boson hopping-pairing matrix in Eq. (2.7) 



acquires an EM Peierls phase: Uf^ — > Uf^ exp[i| f^' ^i^) ■ dr]. Taking a continuum limit, the kinetic part is written 



as 
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where we introduced the boson mass mi, 1/t. The covariant derivative is given by 2> = V + ia^r^ — i^A, where 
we introduced the continuum hmit of the a^j field through af^ = (r^ — rj) ■ a'^ C'"^^^ ). 
The anisotropy term is phenomenologically given in a form 



F^ = ^ Idr 



Cl C3 



(2.17) 



with Cl and C3 being numerical constants of order of unityJll This term describes energy cost associated with small 
fluctuations of the LBC around the SC state {9 = 7r/2). For C3 < ci, the 7-vector prefers to lie in the T1-T2 plane 
(equatorial plane) and the SC state is favored. 
The conventional EM Maxwell term is given by 

The fourth term Fa, the internal gauge field kinetic term, is dynamically induced by integrating out the fermion 
degrees of freedom although we have no such a term initially at the relevant highest energy scales of the fermions 
~ XoJ- We have 

^- = fE E «^('7)n^.(g)«'.('z), (2.19) 



where a = y J A and the fermion polarization bubble originating from the coupling term of the Dirac fermion current 
and gauge field is given by 

n^.(g)=(^..-^)l9|. (2.20) 

We note that this does not take the EM Maxwell form which is proportional to q^, and consequently gives rise to a 
non-local kernel in real space: 

Fa^^ J dr J dr'K{r - r')h{r) ■ h{r'), (2.21) 

where h{r) = V x a^{r) and K{r — r') = er'^'''^'^^'^' ^ Kq with 

nq^l/\ql (2.22) 
instead of Kg = 1 in case of the conventional EM kernel. 

C. London model of a single SU(2) vortex 

In the model of the vortex proposed by Lee and Wen, both a and <j)/2 wind by tt and consequently give an 
appropriate hc/2e vortex for the EM gauge field A{r). This way of winding is specified by 

Va = = ^ (2.23) 
2 2r ^ ' 

which lead to Vtpi = and Viy92 = e^/r, where denotes the azimuthal unit vector in the physical space. That is 
to say, only 62 changes its phase (/?2 by 27r as we go around the vortex, while bi does not. 

The texture of the 7-vector in the SF gauge is indicated in Fig. 1(a). In the SC state outside the core, /j = 
(cos0i, sin0i, 0), while as we approach the core I62I must vanish and the vortex center is represented by li ~ (0,0, 1) 
which is just the SF state. The J^-vector tilts smoothly from the equator to the north pole as the core is approached 



6 



with a length scale denoted by which is identified with the core size. To determine the SU(2) vortex structure, 
we shall use the "London model" prescription of a single vortex in extremely type II BCS superconductor .11 Detailed 
account of the analysis is given in appendix B. 

FIG. 1. (a) The texture of the /-vector in the SU(2) vortex configuration in the SF gauge. At the center of the vortex, 
li points toward the north pole corresponding to the SF state. The shaded circle depicts the vortex core. The local gauge 
transformation gi transforms this configuration to (b) in the d-wave gauge, where the internal phases of the hose condensate 
are gauged away. 



Although quantitative estimation of and is beyond the present simple London model analysis, tc presum- 
ably extends over a fermion coherence length ^ vp/ IS. which may amount to a few lattice scales as suggested 
numericallyH We here just remark that there are two kinds of vortices, because the J-vector can also point toward 
the south pole at the vortex core: 0^ = tt in Eq. ( 2.13| ). This just expresses the state with the staggered flux shifted by 
one unit cell. If the center of the vortex is in the center of the plaquette, the degeneracy between these two kinds of 
vortices is broken by the circulation of the EM superfluid current. This is the situation considered by Wang, Han, and 
Leeil in their numerical local U(l) mean field approach. On the other hand, if the center of the vortex is on a lattice 
site, the degeneracy remains and there is quantum mechanical tunneling between the two states. The tunneling rate 
depends on and is difficult to estimate. However, the dissipation due to quasi-particles may suppress the tunneling 
rate due to the orthogonality catastrophe. Whether the two states are degenerate or not depends on short distance 
physics which is outside the domain of our long wavelength theory. 



III. HOPPING AND PAIRING ORDER PARAMETERS OF THE PHYSICAL ELECTRONS AROUND A 

SINGLE VORTEX 

A. Gauge transformation of the local boson condensate 



Now that the SU(2) vortex model has been established, we shall discuss the effects of the unit cell doubling and 
the phase winding on the hopping and pairing order parameters of the physical electrons around a single vortex. For 
this purpose, it is best to work with the d-wave gauge after making a local gauge transformation by 



The LBC is then transformed to 



gi = exp 





3" 


exp 


I T 

2 



hf - g^hf = e 







(3.1) 



(3.2) 



i.e. the /-vector points toward the north pole everywhere on lattice, as shown in Fig. 1(b). We here consider only 
the case of a single vortex. The great advantage of the d-wave gauge is that the physical electron operator is simply 
written as 



-f 

2 jf'i 



(3.3) 



i.e., the fermions behave as physical electron. After the local gauge transformation to the d-wave gauge, we find 

i/SF ijd _ g.fjSF J 



COS ■ 



sin ■ 



(3.4) 



A,; 



0i + 9i 
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where 



and 



Xij = Acos^^j, Ay = ^sin$y, (3.5) 
$y =$o + (-ir^+-'"%, (3.6) 



4- (3-7) 



As an important consequence of the local gauge transformation, the gauge invariant quantity Vij enters Eq. ( [3.4| ). 
This quantity has a meaning of the fermion "superfluid velocity" associated with the internal gauge field which 



is circulating around the vortex center [see Eq. (B4)]. For example, let us consider Vy{r) along the line, iy = 1/2 



assuming that the vortex center sits at (1/2, 1/2). Fourier transform of Eq. ( BlOj ), gives 



1 1 f°° Jijqi^) 
2r 2 Jo 1 + Aa<7 

where r = [i^, 1/2). In Fig. 2, we show the spatial distribution of Vy{r) by assuming the gauge field penetration depth 
to be Xa = 10 with lattice unit. We see that Vy{r) decays over the length scale Aq, as is naturally expected for the 
superfluid velocity. 

FIG. 2. Spatial distribution of the superfluid velocity Vy{r) associated with the internal gauge field a"'. We assumed the 
gauge field penetration depth to be \a = 10 with lattice unit. The origin should not be taken too literally, since v is defined 
in the continuum limit. 



Let US write Ut^ in the form 



An essential point is that in the d-wave gauge Xij o,nd have the meaning of the hopping and pairing order parameters 
of the physical electron, since the physical electron operator Ci^ is just proportional to the auxiliary fermion operator 



fia [Eq. (3.3)]. Below we discuss the meaning of Xij ^^nd Ay at different limits. 



B. Hopping-pairing order parameters in the vicinity and outside of the vortex core 

1. Vicinity of the vortex center 



First, we consider the vicinity of the vortex center, where 9i ^ 9j ^ and Eq. ( |3.9[ ) becomes 

C74 ^ -At^ cxp [i{-iy^+^y<^>,jT^] , (3.10) 

i.e., 

=v4e^("i)*^^^''*--, (3.11) 
Ay ^ 0. (3.12) 



Eq. (3.12) indicates that the superconducting order parameter is killed at the vortex center. In this region, as is 



directly seen from Eq. (3.11), what is modulated is the phase of the fermion hopping parameter which is just regarded 



as the electron hopping parameter. We see that the sum of the phase around an elementary plaquette yields modulated 
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net flux ±4$o + '/'gauge (j^) with ± signs alternating from plaquette to plaquette. We here introduced a gauge flux 
penetrating an elementary plaquette centered at r 7^ 0: 

0gauge(r) = j( Vxt;.d£^-^/i(r), (3.13) 



where we retained the lattice constant cq. The internal gauge field strength h{r) = V x a'^ is given by Eq. (p3q). 
We see ^gaugeC?^) 4i>o ~ 0{l) [for example, 0(2co) = 0.03 if we take Aq = IOcq]. This situation just indicate the 
fact that the net flux is dominated by the original staggered flux 4$o- Thus, inside the core the staggered phase 
modulation becomes predominant: Ufj ~ — ^r'^exp 1)*"=+-'!' $0^^] • 

The ijfj then breaks not only the translational symmetry (Ufj 7^ t^i+g j+e ' where with fi,!/ — x,y denotes a 
unit vector connecting the neighboring sites ), but also the time reversal symmetry with respect to the local bonds 
{ufj ^ [ufj]*). Although we cannot exphcitly analyze the electronic states inside the core, the time reversal symmetry 
breaking implies that the staggered fermion currents flow on the bonds just as in the case of a uniform SF state. Once 
the bosons are condensed, the currents come up as the staggered orbital currents of the physical hole.0 

In this paper, we are concerned with the possibility of detecting a signature of the unit cell doubling through STM 
measurement. We immediately see that there is no hope in the SF state, because what is staggering in the SF phase 
is the on-bond currents caused by the staggered phase [Eq. (3.11)]. Consequently, the period doubling of the current 
never shows up in the LDOsIi This situation motivates us to look at the region outside the core. 



2. Outside the SF core 



We consider the region outside the SF core. We approximately set 



7r/2 which gives 



Uf 



(3.14) 



I.e. 



Xij — Xij: 



(3.15) 
(3.16) 



Recalling that Xij Aij are interpreted as the hopping and pairing amplitudes of physical electrons, we see that 
the region outside the SF core and inside the gauge field penetration depth, tc ^ f ^ around the vortex, is 
characterized by the staggered modulation of the hopping and pairing amplitudes. Note that the amplitude of Xij ^-nd 



Aij are modulated in a correlated way according to Eq. (3.5) to preserve 



constant. 



(3.17) 



In Fig. 3, we depict the situation given by Eqs. (|3.5[), (3.15), (3.16), and (3.17) 



FIG. 3. Geometric relation of Xij and Aij. The angle ^ij modulates around "l>o in a staggered manner [see Eq. (3.6)] 



ufj breaks the translational symmetry, but does not break the time reversal symmetry with respect to the local 
bonds. Therefore, fjf^ does not cause local fermion current on the bonds [ of course, even in this case, the external 



magnetic field breaks time reversal symmetry and causes globally circulating supercurrent given by Eq. (B6)]. What 
is staggering in this region is not the local current but the local density on the bonds. 

The temporal component of the gauge field, a^^, is determined locally by the LBC. In the uniform case, the saddle 
point is purely imaginary. There, we can regard the LBC as almost uniform in the SC state outside the vortex core. 
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Therefore it may be legitimate to assume a^^ to be uniform and parallel to the LBC (pointing toward the north pole) 
in this region [this assumption is reliable as far as deviation of 9i from ^ 7r/2 is small]. From now on, we set 

= (0,0,ao), (3.18) 

where ao has an order of xJ. 

For the purpose of seeing physical situation, we assume the SF core size to be = 3 and a simple distribution 
of the angle 9 as indicated in Fig. 4(a). In Figs. 4(b) and 4(c), respectively, we show the corresponding spatial 
variation of Xij ^^nd Ay on the link connecting (i^, 0) and (jk, 1) with lattice unit. We also assumed the gauge field 
penetration depth to be Aa ~ 10. Now Ay just represents pairing amplitude of the physical electron and vanishes at 
the vortex center as it should do. As we go away from the core, Xij ^^nd A^ acquire staggered modulation with the 
amplitude becoming smaller, because the superfluid velocity v{r), which is responsible for appearance of the staggered 
modulation, becomes smaller. In fact, the staggered modulation of Xij is just an order of a few percent, while that of 
Aij is rather large. However, as we shall see shortly, period doubling caused by this modulation gives rise to visible 
effects in local density-of-states (LDOS) outside the core. We should also remark that both Xij and A^ contribute 
to the LDOS. The problem now reduces to the more familiar U(l) mean field theory, but with Xij and A^ which 
vary in space. This is precisely the problem treated by Han, Wang, and LeeBiland it is gratifying that they found 
numerically the staggered current around the vortex core as proposed in the SU(2) vortex model.EZI 

FIG. 4. A simple distribution of the 6 indicated in (a) leads to spatial variation of (b) Xij a^nd (c) Aij on the link connecting 
{ix, 0) and (ix, 1) as indicated in the inset. The gauge field penetration depth is assumed to be Aa = 10. 

In Fig. 5, we schematically show the modulation pattern of Xij outside the core. The staggered modulation becomes 
most conspicuous when scanned along the straight line ix — 1/2 or iy = 1/2, provided that the vortex center sits 
at (1/2, 1/2), because on these bonds the circulating 7;(r)-ficld becomes parallel to the bond directions. Apparently, 
the bond modulation pattern reminds us of the spin-Peierls states. However, this is not the case, since the MF 
expectation value of spin-exchange energy on the bonds is given by {Si ■ Sj) — ~J{Xij + ^fj) = const, and therefore 
the spin-Peierls order parameter becomes {Si ■ Sj+ej, ~ Si ■ Si-e^) = with being a unit vector connecting the 
neighboring sites. 

FIG. 5. Schematic drawing of the amplitude modulation pattern of the hopping parameter Xij outside the SF core. Solid 
and dotted bonds indicate enhanced and reduced amplitudes, respectively, where thickness of the bonds qualitatively represents 
magnitude of the modulation. Circulation of the the fermion "superfluid velocity", v{r), associated with the internal gauge 
field is indicated by the arrows. The staggered modulation becomes most conspicuous when scanned along the lines ~ 1/2 
or iy = 1/2, provided that the vortex center sits at (1/2, 1/2). Note that the boundary of the SF core region, inside which the 
staggered orbital currents flow, should not be taken literally. In reality, there is a crossover region around r ~ where the 
staggered current and the staggered amplitude modulation coexists. 

As we approach the core from the outside, the /-vector in the SF gauge gradually rises off from the equatorial 
plane [see Fig.l (a)]. This may give rise to a crossover region characterized by coexistence of the amplitude and 
phase modulation, where 0-dependence of Ufj becomes significant. It is expected that the staggered current begins to 
appear around r ic and its strength becomes stronger as we approach the immediate center of the vortex. We give 
a schematic drawing of this circumstance in Fig. 5. To study the effects of ^-dependent Ufj is, however, beyond the 
scope of the present paper and we concentrate on the region £c r Xa- We should also remark that when the angle 
6 deviates from di ~ 6j = tt/2 a.s we approach the core, the direction begins to slightly deviate from the north 
pole, since is no longer parallel to the /-vector due to small anisotropy. In the next section, we shall compute the 



LDOS in the SC state outside the core by setting 9i = 6j = n/2. Then, a^^ is given by Eq. ( 3.18 ) and exactly parallel 
to the /-vector pointing toward the north pole. We expect our results to be qualitatively valid even for r ~ as long 
as we avoid the inside of the core. 
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IV. LDOS OUTSIDE THE CORE 



As we saw in the previous section, the staggered modulation of the hopping amphtude Xij and pairing amphtude 
Aij becomes predominant over the region ic ^ f ^ ^a- The presence of staggered modulation suggests that this may 
be the best place to look for unit cell doubling effect. In this section, we consider the LDOS in this region. 



A. Formulation of LDOS 



The local density of states at an arbitrary point r on lattice is defined by 



N{r, Lo) = --ImgP'-y^r, r; tuj) U^=^+,s , (4.1) 

IT 

where the propagator for the physical electron is introduced by G^^^^ir, r', iuj) = — dre"^'^ < TrCa-{r, r)cj.(r') > . 
To model the tunneling current we assume that the electron tunnel from the tip located at r to a linear combination 
of Wannier orbitals centered at lattice sites, i.e., the physical electron operator at r, Ca{r), is related to Cia as 

c-W =II"»We''^^'''''''''^^'''^c,,> (4-2) 

i 

where the EM gauge potential A gives rise to the EM Peierls phase. The envelope function ai{r) may be simulated 
by ai{r) — e^l''^''''/^ in the bond direction (the Cu-O-Cu bond). The length scale ^ can reasonably be set equal 
to ^ — 1/2 with lattice scale, corresponding to the Cu-0 separation. Since the effects of the EM gauge fields are 
negligibly small in strength as compared with the internal gauge potential, from now on, we ignore the EM Peierls 



phase and examine the effects of the staggered hopping and pairing amplitudes on the LDOS. Noting Eq. (3.3) in the 



d-wave gauge, Eq. (4.1) is written as 



N{r,uj) = Im^at{r)aj{r)[gfj{iuj)]ii\i^^^+is. (4.3) 

The subscript 11 means 11-component of the lattice fermion propagator of a 2 x 2 matrix form, Gfj{T) = 

We here give an intuitive demonstration that the LDOS exhibits conspicuous staggered pattern only when measured 
on the bonds. More quantitative discussion will be given in the following subsections. For example, we pick up the 
sites 1,2, ...,6 indicated in Fig. 6 and consider the midpoints on the bond, Bi, B2, and the plaquette centers, Ci, C2. 
The LDOS at Ci and C2 come from X)?; 2,4,5 ^ij"' Z^i j=2, 3,5,6 ' respectively. We see, however, ^ QIq 
because the bonds 12 and 56 are almost equivalent except the effects of negligibly small dependence of the v-field on 
the spatial position r over the lattice scales. Similarly, Q^^ ^ C/fg and Q^/^ ~ C/fg. Therefore, N{rci,uj) ~ A^(rc2,'^)- 
Similarly, the LDOS at the lattice sites is almost uniform. On the other hand, the LDOS at Bi and B2 come from 
J2i j=i 2 Sfj' ^^'^ X^i j=2 3 Sfj ' respectively. Here, and are clearly inequivalent because they connect the bonds 
with alternating hopping-pairing amplitudes. 

FIG. 6. Points on lattice where we consider the LDOS. We have four symmetrically distinct points: the plaquette-center 
(x), site-top (•), and bond-center (o). The site-top and bond-center points correspond to the Cu and O sites, respectively. 

To compute Gfj (ioj) in the SC state outside the core, we shall use the following two approaches which may be com- 
plementary to each other: (I) a perturbative analysis using the gradient expansion , and (II) an exact diagonalization 
using the "uniform v "approximation . In the former approach, we can take account of the circulating configuration 
of the t;(r)-field, while in the latter approach, instead, we can obtain a nonperturbative aspect of the problem. 
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B. Perturbative analysis using the gradient expansion 



First, we expand (3.5) with respect to Vij up to the first order as 

X^J ~ Xo - {-iy^+'-Aov,„ (4.4) 

^ Ao + i-ly^+^yXov^J, (4.5) 

which give C7^ = f/^ + 6Uf^ with 

SU,, = (-1)'^+^«[Aot3 + (^lyy+^yxoAv.,. (4.6) 

Then, we treat the term 

•^^^ = i E 45U.3^^. (4.7) 

<ij> 



as perturbation with respect to where aoi is given by Eq (3.18). The free propagation is governed by Uf^ and the 
corresponding propagator becomes 

eo"(fe,-) = -^ + -^, (4.8) 

where the generahzed coherence factors are introduced by ?7fc = ^[l + {ikT^ + 'HkT^) I Ek\, and 14 = ^[1 — 
[ikT^ + ilkT"^)/ Ek]- The one-particle spectrum is given by 



Ek^^ll+ril (4.9) 

with 7fc = —Jxo[coskx + cosfcj, + ^2 cos fc^; cos fcj^ + t^{cos2kx + cos2fcj^)] + oq, and r/^ = +JAo(cosfc2. — cosfcj^). 
We took account of the second and third nearest neighbor hopping of the fermions to reproduce the real 
band structure. In general, the d-wave nodes shift from (±7r/2, ±7r/2). In the case of t2 = t-^ — 0, the 
nodes are located at (± cos~^[ao/2 Jxo]i i cos~^ [ao/2 Jxo])- For ^2 7^ and ^3 7^ 0, the nodes are located at 
(± cos~i[/(t2 J3, ao)], ± cos~^ [/(f2 J3, ao)]), where 

/(i2, ta, ao) = ^- . (4.10) 

413 + 12 

We see that as far as ^3 7^ the nodes shift from (±7r/2, ±7r/2) even if oq = 0. Furthermore we note that location of 
the nodes is independent of the gap magnitude A. 

Since the perturbation term causes period doubling, it is convenient to introduce the fermion operators on two 
sublattices: 

V', = ^ E (^-^ ± ^^+q) ' (4.11) 

where Q = (tt, tt) and k g RZ means k runs over the reduced Brillouin zone + |fcj,| < tt. We dropped the spin 
indices. The + and — signs are for the cases where i belongs to the A [r^ = [i^, iy) = (even, even) or (odd, odd)] and B 
[{ix,iy) = (even, odd) or (odd, even)] sublattice sites, respectively. Then, as derived in appendix C, the perturbation 
term is written in momentum space as 

^-Y. Et^I+f +Q.^'^(9)^'^-§- + H.C.], (4.12) 

fcSRZ q,cr 
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where Ck{q) = AoC+{q)T^ + xoC^ {q)T\ with 

Cfc (9) = ttJ ^ M^Ll bysmfc^ ig^sin/cj,] . (4.13) 

\q\ L + Aa\q\ 

The momentum transfer q should be small because we retained only slowly varying v-field. The perturbation processes 
cause the unit cell doubling and scatter the electron with k in the reduced zone to A; + Q in the second zone, and 
consequently the mirror image of the reduced zone is formed in the second zone, as indicated in Fig. 7(a). 



FIG. 7. (a) The perturbation processes given by Eq. (4.7) connect the electron with k in the reduced zone [inner square] 
to fe + Q in the second zone [shaded region], and consequently the mirror image of the reduced zone is formed in the second 
zone. The rf-wave nodes inside the reduced zone and their mirror images are also indicated, (b) The scattering processes along 
the (0,0) (7r,7r) direction whose matrix elements give the coherence factors and A^^. At the energy w* , the Dirac cones 
around the d-wave nodes touch the reduced zone boundary and the resonance occurs, (c) The level crossing at the reduced zone 
boundary would be lifted and eventually the period doubling would cause gap opening if we would go beyond the perturbative 
scheme. Note that situations in (b) and (c) correspond to the case of a simple band structure without the next (12) and second 
nearest (ta) fermion hopping. 

Now, we consider the four distinct points on lattice indicated in Fig. 6: (a) the center of the plaquette (plaquette- 
center), (b) the top of the sites (site-top), and (c) the center of the bonds (bond-center). The site-top and the 
bond-center points correspond to the Cu and the O sites, respectively, on the Cu02 plane. All the detail of derivation 
of the LDOS is left to appendix C. In any case, the LDOS is written in a form 

N{r,u)/xa^ ^ Na{Lo) ± 6N{r,uj), (4.14) 

where -|- and — signs alternate from plaquette to plaquette, site to site, or bond to bond for the cases (a), (b) and (c), 
respectively, and a represents magnitude of the envelope function from the nearest site. The uniform counterpart is 
given in a form 

Noicu) = --lmy^^Uk)[g^{k,k,^Lo)]n\^^^u+^s, (4.15) 

k 

where we introduced the generalized propagator Q^{k,k',iu;) = r ''^''^i^ (*'^)- matrix elements 

Mo(fc) distinguishes different symmetries associated with each point and are given by 

^^plaqucttc^^^ = C0S2 COS^ ^, (4.16) 

M^^^k) = 1, (4.17) 
Mo^°"'i(fc) =cos2y, (4.18) 

for the cases (a), (b), and (c), respectively. The perturbation processes do not affect the uniform counterpart within 



the Born approximation, and thus, in Eq. (4.15), we obtain 



lm[g {k,k,iuj)]ii\i^^^+is ^lin[gQ{k,iuj)]ii\i^^^+is (4.19) 

which just reproduces the LDOS prohlc in the uniform d-wave SC state except overall reduction due to the matrix 
element Mo{k). 

The staggered counterpart is given in a form 

5N{r,Lu)^-hm ^ ^ Af (fe, 9; r)[5F(fc + | + Q, fc - |, ic^) + ^^(fc _ |, fc + | + Q, ic^)]ii (4.20) 

q-^ small fc^RZ 
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The matrix elements M{k, q; r) associated with each point are given by 

^^piaqucttc(^^ q; r) ^ cos{q ■ r) sin ^^-^ sin ^^^^ cos cos (4.21) 

Af"*'=(fc,q;r) = cos(q-r), (4.22) 

M'^°"'^(fc, q; r) = sin(g • r) sin ^ ^ cos ~ ^ , (4.23) 

where r denotes the plaquette-center, site-top, and bond-center points, respectively. Now we need to compute 
G^{k + ^ + Q,k — ^,iuj) and Q^{k — ^,k + ^ + Q,iuj). The detail of computation is presented in appendix C. 
We obtain 



[^F(fe+| + Q,fc- |,iLj) + gF(fc- |,fe + |+Q,iw)]n U^^^+u 



+ |5(c., -E^+^/2+Q,+Ek-^^2) [AoC+{q)N^- + XoCfe (g)7Vr] , (4-24) 



where S{lu, x, y) = [S{lu — x) — S{ll! — y)]/{x — y). The first equality in Eq. (4.24) just represents the inversion and time 
reversal symmetries of the propagator: t/^(fc + f + Q, fe — f , ^i^) = ^^(fc — f , + f + Q, iw). The coherence factors 
are given by 

= 1 + 7+7- - V+Tl- ± 7+ ± 7-, 
^fc"^ = ±^7+ ± ?7- + l+V- + '7-H7-, 
^i^"^ = 1 - 7+7- + -n+V- ± 7+ T 7-, 

^fc = ±'7+ TV-- 1+V- - V+1- , (4.25) 

where 7+ = 7fc+q/2+Q/-E^fc+q/2+Q, 7- = lk-q/2/ Ek-q/2, V+ = Vk+q/2+Q/Ek+q/2+Q, and 77- = ■qk-q/2/Ek-q/2- 

To proceed with further analytical computation, we note that main contribution of k integral comes from regions 
near the nodes in the vicinity of (±7r/2, ±7r/2), while q is small. Thus, it is legitimate to ignore q with respect to k 
in Mo{k), M{k,q]r), and G^ik + ^ + Q,k~ ^,iu>), while we must retain q in C^{q). This approximation amounts 
to ignoring r-dependence of t;(r)-ficld over the lattice scales, and retaining only fermion fluctuations. On the other 
hand, retaining q-dependence of C^(<7) amounts to taking account of long distance decay of the t;-field. Under this 
approximation, Eqs. (4.21), (4.22), and ( 1.23| ) are simply reduced to 



^piaquctto^^^ q;r) = ^ cos(q • r) sink^ sinfcy, (4.26) 
A'r"°(fc, q; r) = cos(q • r), (4.27) 
M^°'''^{k, q;r)^^ sm{q ■ r) sin k^. (4.28) 

By noting the antisymmetry relation C^{q) = —C^{—q), we immediately see that SN{r,uj) vanishes at the plaquette- 
center and site-top points while it remains finite at the bond-center points. Thus, we confirm that the staggered 
counterpart of the LDOS appears only when measured on the bonds. Even in the cases of the plaquette-center and the 
site-top, dN{r,uj) becomes finite if we retain q with respect to k, i.e., take account of negligibly small dependence of 
the v-field on the spatial position r over the lattice scales. However, this effect is still invisibly small as compared 
with the case of the bond-center. This result is fully consistent with an intuitive discussion given in Sec. IV A. 
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From now on, we concentrate on the bond-center points: the midpoint of the bond connecting i and i + where 
fi — X or y. Taking account of the envelope function, the magnitude of the LDOS may be reduced by a factor 



e ^ ~ 0.1 as compared with the uniform counterpart of the LDOS at the site-top. Using Eqs. (4.24) and (4.28), the 
q-integration in Eq. ( 4.20| ) can be performed to yield 



SN{r,u;) = LJ}^y^(^r) ^ sin^ A:^ (4.29) 

+N+S{uj; Ek, ~Ek+Q) + N^S{uj; ~Ek, Ek+q)] , 

where the coherence factors = AqL'^'^ + XoLk"^ , and = AqN^"^ + Xo^fc^ represent the matrix element 
associated with the scattering processes indicated in Fig. 7(b). As has already been mentioned, the best paths to 
detect the staggered modulation of the LDOS are the lines Za; = 1/2 or ij, = 1/2 provided that the vortex center sits 
at (1/2,1/2), because in this case we can go through the bonds whose directions are parallel to the circulating 
tj(r)-field [see Fig. 5]. 

1 . The case of t2 — — 



First, we consider a toy band structure with ^2 = ^3 = in Eq. (4.9), because this simple case provides us with a 
clear view on the period doubling effects. In Fig. 8(a), we show the profile of N{r,uj)/xa'^ at the four bond-center 
points, A{ix, 1/2), B{ix + 1/2,0), C{ix + 1, 1/2), and D{ix + 1/2, 1) with ^ 5 [see the inset of Fig. 8(a)]. From now 
on, we fix the parameters oq = 0.05xo>^: ^o/xo — 0.2, and assume the gauge field penetration depth to be Aa = 10. 
This choice of uq and Ag is reasonable in the underdoped regime.B Note that at B and D, 6N{r,uj) almost vanishes 
and the LDOS is just given by No{lu), because v{r) becomes almost perpendicular to these bond directions. The 
modulation pattern at the other points can be read off from Fig. 5. 

FIG. 8. (a) LDOS profile in the case of t2 = ts = 0, obtained by the perturbative analysis at the points A, B, C, and D 
indicated in the inset. The LDOS at B and D are just Nq{uj). The peaks at to = u/xoJ ~ ±0.38 are associated with the d-wave 
superconducting gap. The additional peaks at Cj = ±0.41 are associated with the van-Hove singularity located at (0, ±7r) and 
(±7r,0). The staggered structure around u = 0.05 comes from resonant scattering between the fermions with k and k + Q, 
caused by the period doubling, (b) The one-particle energy contour around the d-wave node, (c) The energy contours Ek ~ u> 
and Ek+Q ~ touch at d; = ±0.05. 

We see that inside the overall V-shaped profile with the sharp peaks at a) = uj/xqJ ~ ±0.38 associated with the 
d-wave superconducting gap, there appears additional peak and dip structure at site C and A, respectively, around 
Lj = ±0.05. From now on, we refer to this structure as "staggered peak-dip (SPD)" structure, since the peak and dip 
alternate from bond to bond in a staggered manner. The additional peaks at d) = ±0.41 come from the van-Hove 
singularity located at (0, ±7r) and (±7r, 0) points. The low energy dispersion gives elliptic contours around the d-wave 
node as indicated in Fig. 8(b) which touch the reduced zone boundary at (7r/2, 7r/2) as the energy increases. The 
specific structure around ui = 0.05 comes from resonant scattering between the fermions with k and k + Q. As cj 
increases from zero, the energy contours Ek — oj and Ek+Q = lo touch at (7r/2, 7r/2) on the reduced zone boundary at 
Cj — ±0.05 as indicated in Fig. 8(c) [see also Fig. 7(b)] and resonance occurs. We note that the modulated structure 
inside the V-shaped profile is predominant on the particle side (w > 0). This asymmetry is due to the matrix element 
effect: vanishes at (7r/2,7r/2). 

In any case of this toy band structure, it may be totally hopeless to experimentally detect such tiny structures as 
indicated in Fig. 8(a). We see in the following that the realistic band structure of BSCCO drastically changes this 
situation. 
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2. The case of real hand structure 

Next, we take account of t2 — —0.550 and ^3 = 0.087 to reproduce the real band structure of BSCCO measured by 
angle-resolved photo emission spectroscopy.il In Fig. 9(a), we show the profile of N{r,uj)/xa^ at the same points as 
in Fig. 8(a). In this case, inside the overall V-shaped profile with the sharp peaks at J) = oj/xqJ ~ ±0.323 associated 
with the d-wave superconducting gap, there appears prominent SPD structure around uj = ±0.179 and uj — 0.224. In 
Fig. 9(b), we show the same profile as in Fig. 9(a) over an wider energy window. The peaks ai uj = ±0.79 are ascribed 
to the van-Hove singularity at (0,±7r) and (±7r,0) points.Ei 

The SPD structure inside the V-shaped profile again comes from resonant scattering between the fermions with k 
and k + Q. As seen in Fig. 9(c), the low energy elliptic contours in the case without t2 and ^3 [Fig. 8(b)] bend around 
the d-wave nodes [bending of the Dirac cone]. Consequently, as lo increases from zero, the energy contours = lo 
and Ek+Q = lu first touch on the reduced zone boundary at lu — ±0.179 as indicated in Fig. 9(d) and resonance 
occurs. Then, at cj — ±0.22 they touch again at (7r/2, n/2) as indicated in Fig. 9(e) and the second resonance occurs. 
The reason why the second resonance comes up only in the electron {cu > 0) side is again ascribed to the matrix 
element effect as in the case of t2 — — 0. We can say that due to the real band structure [bending of the Dirac cones 
around the d-wave nodes] the staggered structure in the LDOS profile becomes far more prominent as compared with 
the case of t2 = t^ = 0. 

We see that the SPD structure due to the period doubling occurs only inside the V-shaped profile [see Fig. 9(b)]. In 
fact, the energy scale at which the SPD structure appears depends on the band structure parameters [ao oc x, Aq/xo, 
t2, and ts]. For a reasonable choice of the parameters in the underdoped regime, however, the resonance always occur 
at the energy scales below that of the superconducting gap, i.e., the SPD structure always appear inside the V-shaped 
profile. 

FIG. 9. (a) LDOS profile in the case of the real band structure of BSCCO, obtained by the perturbative analysis at the 
points A, B, C, and D indicated in the inset. The LDOS at B and D are just No{u)). The staggered structure around ti — ±0.179 
and Cj = 0.224 comes from resonant scattering between the fermions with k and k + Q, caused by the period doubling. The 
small wiggles outside the V-shaped profile come from numerical fluctuations, (b) The profile over a wider energy window than 
that of (a). The peaks at = ±0.79 are ascribed to the van-Hove singularity at (0, ±7r) and (±7r, 0) points, (c) The one-particle 
energy contour around the d-wave node. The energy contours Ek — ui and E^+q ~ u) touch at d; = ±0.179 and ±0.224 as 
indicated in (d) and (e), respectively. 

To see a qualitative feature of the doping dependence, in Fig. 10(a), we show the LDOS profile for cq = 0.03xo>^ 
and Ao/xo = 0.35, corresponding to the case of a lower doping as compared with the case of oo — 0.05xo'^ and 
^o/xo — 0.2. We see that the SPD structure remains robust, although the resonance occurs only once at w = 0.2. 
Smearing out of the second resonance is due to change of the geometry of the Dirac cone around the c?-wave nodes. 
The shape of the low-energy contours change upon changing A as clearly seen by comparing Fig. 10(b) with Fig. 9(c). 
The contours E^ = uj and E^+q = uj touch on the reduced zone boundary only at uj — ±0.2 [Fig. 10(c)]. As already 
mentioned, however, location of the d-wave nodes is independent of Aq and always shift from (±7r/2, ±7r/2) for finite 
is, i.e, the resonance at (±7r/2, ±7r/2) occurs at the enregy 

(i* =±(2i3±ao/xoJ). (4.30) 

In this respect, the next nearest neighbor hopping t^, plays a crucial role to push the energy scales of the SPD structure 
toward visibly finite energy scales. 

FIG. 10. (a) LDOS profile for ao = O.OSxoJ and Aq/xo = 0.35, corresponding to the lower doping as compared with 
the case of ao = O.ObxoJ- (b) The one-particle energy contour around the d-wave node. The energy contours Ek = uj and 
Ek+Q = UJ touch at = ±0.2 as indicated in (c). 
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In the perturbative picture presented here, the period-doubled perturbation processes form the "mirror image" of 
the energy bands with respect to the reduced zone boundary [Fig. 7(a)]. The energy levels Ek and its mirror image 
Ek+Q cross on the zone boundary \kx \ + \ ky\ — it, which causes the resonant scattering at the corresponding energy 
uj* [Fig. 7(b)]. It is naturally expected that if we go beyond the perturbative scheme the level crossing would be 
lifted and eventually the period doubling may cause gap opening in the fcrmion excitation spectrum as indicated in 
Fig. 7(c). This point is confirmed through the exact diagonalization under uniform v approximation as shown below. 



C. Exact diagonalization using the uniform v approximation 



Next we consider the case of uniform u-field: = {v^xi voy) which may locally capture the effects of the circulating 
v{r). From Eqs. (3^), ( 3.15| ), and ( 3.16 ), we see that the uniform Vq yields 

X^J = X^] = ^cos[$o + (-1)''"'+^''' (r, - r j) ■ vo], (4.31) 
A„- = Ay = ^sin[$o + i-iy-+'«in - r,) ■ Vo]. (4.32) 

An advantage of the uniform v approximation is that we can exactly diagonalize the corresponding fermion Haniilto- 
nian, which can be written as 



^0 = ^ E *LTfe*fe,, (4.33) 

feGRZ 

where {^kaV = ((V'fccr)^, {ipk+QaV)- The 4x4 matrix Tfc is given by 



Vk + akT'^ iWk 

iWk+Q Vk+Q + ak+QT 



(4.34) 



where Vk = -XqJ%t^ + AoJflkT^, and Wk = AoJAfcr^ + Xo^MfcT\ with % = cos uqx cos fc^; + cos VQy cos ky,rjk = 
cos vqx cos kx — cos VQy cos ky, Xk = sin vqx sin kx + sin VQy sin ky, flk — sin vqx sin kx — sin v^y sin ky . Noting the fact that 
the field v does not modulate the hopping amplitude between the same sublattice sites, we take account of the hopping 
parameters t2 and by introducing 

dk — 0.0 ~ ^2 cos kx cos ky — i3{cos2kx + cos2ky). (4.35) 

The one-particle propagator in a 4 x 4 matrix form is given by 

G^{k,iuj) ^[iujl-TkV^, (4.36) 

where 1 denotes a 4 x 4 unit matrix. As was inferred from the perturbative analysis, the unit cell doubling brings 
about the one-particle spectrum split into two branches in the reduced zone, ±i?^ and i-E;^, where 

Ei = [al + A'fc + Mfc ± 2{a2 (7^ + ) + (ry^A^ + jkf^kfV^V^ (4-37) 

with k e RZ. To compute the LDOS, we need 11, 33, 13, and 31 components of G^{k,iijj) which are explicitly given 
in appendix D. 

As in the perturbative analysis, we consider the LDOS at four distinct points on lattice: (a) the plaquette-center, 
(b) the site-top, and (c) the bond-center points. Repeating an analysis similar to that in appendix C, we obtain the 
LDOS in a form 

N{uj)/xa^ = iVo(w) ± 6N{uj), (4.38) 
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where + and — signs alternated plaquette to plaquette, site to site, or bond to bond for the cases (a), (b), and (c), 
respectively. The uniform counterpart is given exactly the same form as in the case of the perturbative analysis: 

^o(^) = --Im VAfo(fc)[GF(fc,iw)]n |„+,5, (4.39) 

TT ' * 



where Mo(fc) are given by Eqs. (4.16)-(4.18). We here used the relation, [G^(fe, ici;)]33 — [G^(fc + Q, icj)]ii, which is 
explicitly shown in appendix D. 

The staggered counterparts at the plaquette-center and the site-top points are given in a form 

(5iV(w) = -ilm ^-^W ([G^(fe,«^)]i3 + [G(fc,zc^)]3i) (4.40) 



TT 

fceRZ 



where the matrix elements M{k) are given by 

^plaquette (j^) = sink^ sinfcy, (4.41) 
iVr'*'=(fc) = 1. (4.42) 

As shown in appendix D, we have the following relation: [G^ {k,iuj)]i3 — — [G^(fc, zcj)]3i = [G^ {k,iuj)]'^i. Therefore, 
6N(lu) exactly vanishes at the plaquette-center and the site-top points. 
On the other hand, at the bond-center points we obtain 

SNiio) = -hm sinfc^(i[GP(fc,iw)]i3-i[G(fc,ic^)]3i) (4.43) 
fceRZ 

which remains finite, where we considered the bond in the e^-direction. Thus, just as in the perturbative analysis, the 
LDOS exhibits staggered pattern only when measured on the bonds. Using an explicit form of G^ given in appendix 
D, we obtain 

Noico) = ^cos^ ^Ukico) [6iiu - E^) + 5{u: + E^) - 5{u - E+) - 5{u, + E+)\ , (4.44) 
fe 

5N{u)= Y KUk{cj)sinkf,[S{u;-E^) + S{u; + E^)-S{uj-E+)-S{u; + E+)], (4.45) 
fceRZ 

where Uk{uj) = [w^ -I- LoAk ~ Bk - Ck/uj]/2Pk, and Ukito) ^ -[lu + 2au + Ckl'AI'^Pk-, with Pfc, A^, ... being given in 
appendix D. 



1. The case of t2 ~ — 



First, we consider again a toy band structure with <2 = ^3 = in Eq. (49). In Fig. 11(a), we show the profile of 
No{u!) and Nq{u!) ± SN{llj) for Vq = (0,0.1), of which direction and strength locally simulate v{r) around the points 
B, D, and A,C in the inset of Fig. 8(a), respectively. We used the same parameter set as in the case of Fig. 8(a). 
In Fig. 11(b) is indicated the energy contour of the lower band E^ with the corresponding band structure of E^ 
being shown in Fig. 11(c). The uniform Vq field breaks the original 4- fold symmetry and the d-wave nodes are located 
slightly off the F— M line. The van-Hove singularity on the Y-F line is caused solely by the superconducting gap and 
gives peaks at a) = ±0.38. 

FIG. 11. (a) Profile of iVo(tj), and No{iu) ± 5N{iu) for the uniform field vo = (0,0.1) in the case of t2 = U = 0. In the 
inset is shown fine structure of Nq{ijj) around Co ~ 0.1, detected with higher numerical resolution, (b) The energy contour of 
the lower band E~ and (c) the dispersion of with vq = (0,0.1) along the path r(0, 0) M(7r/2, 7r/2) Y(0,7r) V. 
Fine band splitting on the reduced zone boundary are magnified in the inset, (d) Profile of No{u}), and No{ijj) ± 5N{uj) for the 
uniform field vq — (— 0.1/\/2, 0.1/-\/2). In the inset is shown fine structure of iVo(aj) around a) ~ 0.05, detected with higher 
numerical resolution, (e) The energy contour of the lower band Ej^ and (f) the dispersion of E^ with i;o = (— 0.1/\/2, 0.1/\/2). 
Fine band splittings on the reduced zone boundary are magnified in the inset. 
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As expected from the perturbativc analysis, the unit cell doubling causes the gap opening on the reduced zone 
boundary between ui = 0.10 and oj = 0.12. The van-Hove singularities associated with this gap structure gives rise to 
the specific structure in the LDOS profile. The corresponding fine structure in No{ui) could be detected with much 
higher numerical resolution[720x720 meshes of the Brillouin zone], as shown in the inset of Fig. 11 (a). The van- Hove 
singularity at Y(0,7r) point is intrinsic to the normal state dispersion 7/5 and give peaks at a; = ±0.41, just as in the 
case of Fig. 8(a). 

We see that the staggered modulation profile shown in Fig. 11(a), iVo(w) ± SN{uj), is in remarkable agreement 
with Fig. 8(a) obtained by the perturbativc analysis. However, a striking difference is that the dip structure around 
OJ = 0.105 is now intrinsic to the modified band structure with van-Hove singularities associated with the gap opening 
on the reduced zone boundary and appears even in the uniform counterpart No{u!). 

In Fig. 11(d), we show the LDOS profile for vq = (— 0.1/\/2,0.1/\/2), where strength of vo is the same as in the 
case of Fig. 8(a), but its direction locally simulates v{r) in the 45 degree direction in Fig. 5. In Fig. 11(e) is indicated 
the energy contour of the lower band with the corresponding band structure of being shown in Fig. 11(f). 
The uniform Vq = (— 0.1/\/2j0.1/V2) field breaks the original 4- fold symmetry in a way different from the case of 
Vo = (0,0.1). Consequently, the energy scales of the van-Hove singularities responsible for the dip structure move 
downward. The qualitative feature of the profile, however, does not change much for small magnitude of vo considered 
here. 

2. The case of real band structure 

Next, we turn to the real band structure of BSCCO. In Fig. 12(a), we show the profile of iVo(cj) and No{^) i SN(uj) 
for Vo = (0,0.1), of which direction and strength locally simulate v{r) around the points B, D, and A,C in the inset 
of Fig. 9(a), respectively. We used the same parameter set as in the case of Fig. 9(a). In Fig. 12(b) is indicated the 
energy contour of the lower band E^ with the corresponding band structure of E^ being shown in Fig. 12(c). The 
van-Hove singularity on the Y-F line is caused solely by the superconducting gap and gives peaks at a) = ±0.323. 
The uniform Vq field breaks the original 4-fold symmetry and the d-wave nodes are located slightly off the F— M line. 
The unit cell doubling causes the gap opening on the reduced zone boundary between uj — 0.229 and uj — 0.265, 
corresponding to the second resonance in the perturbativc analysis [the second touch of the energy contour indicated 
in Fig. 9 (d)]. In this case, due to the presence of t2 and additional van-Hove singularity occurs at a) = 0.186 and 
u) = 0.216, corresponding to the first resonance in the perturbative analysis [the first touch of the energy contour 
indicated in Fig. 9 (c)]. The corresponding fine structure in No{^) could be detected with much higher mmicrical 
resolution[720x 720 meshes of the Brillouin zone], as shown in the inset of Fig. 12 (a). As was mentioned in the 
perturbative analysis, the van-Hove singularity at Y(0, tt) point is pushed upward as compared with the case of 
^2 = ^3 = and lies at the energy Co = ±0.79 outside the energy window of Fig. 12(a) just as in Fig. 9(b). 

FIG. 12. (a) Profile of No{iL>), and No{u}) ± SN{u>) for the uniform field Do = (0,0.1) in tlie case of tlie real band 
structure of BSCCO. In the inset is shown fine structure of No{u>) around lj ^ 0.2, detected with higher numerical res- 
olution, (b) The energy contour of the lower band Ej^ and (c) the dispersion of with vo = (0,0.1) along the path 
r(0, 0) — > M(7r/2, 7r/2) Y(0, tt) F. Fine band splitting on the reduced zone boundary are magnified in the inset, (d) 
Profile of No{(^), and Noi^o) ± SN{uj) for the uniform field vo = (— 0.1/\/2, 0.1/\/2). In the inset is shown fine structure of 
Nq{ijj) around uj ~ 0.2, detected with higher numerical resolution, (e) The energy contour of the lower band and (f) the 
dispersion of E^ with vo = (-0.1/^,0.1/^2). 

We see again that the profile in Fig. 12(a), in particular the SPD structure inside the V-shaped profile, is in 

remarkable agreement with Fig. 9(a) obtained by the perturbative analysis. However, as in the case of ^2 = ^3 = a 
striking difference is that the SPD structure is now intrinsic to the modified band structure with van- Hove singularities 
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associated with the gap opening on the reduced zone boundary and appears even in the uniform counterpart Nq{uj). 
This suggests that in reahty the SPD structure may be detected not only on the bonds but also at sites. 

In Fig. 12(d), we show the LDOS profile for Vq = (— 0.1/\/2, 0.1/-\/2), where strength of Vq is the same as in the 
case of Fig. 9(a), but its direction locally simulates v{r) in the 45 degree direction in Fig. 5. In Fig. 12(e) is indicated 
the energy contour of the lower band E^^ with the corresponding band structure of being shown in Fig. 12(f). 
Qualitative feature of the LDOS profile given in Figs. 12(a) and 12(d) are quite similar. Thus, we may say that the 
SPD structure is robust and detectable in all the directions around the vortex center. 

We note that in both perturbative and exact analysis, the SPD structure in the LDOS is predominant on the 
particle side (w > 0) . We can understand this asymmetry by first turning off the superconductivity and consider the 
effect of unit cell doubling. Since we are doping with holes, the gaps being opened by the unit cell doubling are on 
the empty side on the Fermi surface. Matrix element effect preserves this particle-hole asymmetry even after we turn 
on the superconductivity. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have concentrated on how to detect a signature of the unit cell doubling originated from the SF 
state through STM measurement. Although the signature of the SF state appears only dynamically in a uniform SC 
state, a topological defect (vortex) stabilizes static texture of the boson condensate and the spatial component of the 
massless internal gauge field . We determined the texture associated with a single vortex based on a simple London 
model. A half flux quantum of the EM gauge field. A, penetrates over a huge region r < A^, as compared with a 
half flux quantum of the internal gauge field, a^, which penetrates over a region r ^ Xa- Although the fermions do 
not couple to the EM gauge field, they still see the internal gauge flux tube associated with a^. By this reason, the 
topological texture shows up in the hopping (xij) a-nd pairing (A^) order parameters of the physical electrons and 
gives rise to the staggered modulation of Xij ^^'^ ^ij through the gauge invariant "superfluid velocity"!) associated 



with [see Eqs. (gj) and ( p. 61 )] 



The most important formula in this paper is Eq. (3^) which directly tells us that whereas the center in the vortex 
core is a SF state, as one moves away from the core center, a correlated staggered modulation of the hopping amplitude 
Xij and pairing amplitude of the physical electrons becomes predominant over the region ic ^ r Xa- Combining 
the results obtained through the gradient expansion and the uniform v approximation, we concluded that the signature 
of the unit cell doubling may be most prominently detected through the staggered peak-dip (SPD) structure inside the 
V-shaped profile measured on the bonds. The real band structure of BSCCO, in particular the next nearest neighbor 
hopping ts, plays a crucial role to push the energy scales of the SPD structure toward visibly finite energy scales. 
The structure directly originates from the unit cell doubling, which is stabilized by the topological texture (phase 
winding) under the external magnetic field. In this respect, our effects have little to do with the d-wave symmetry of 
the superconducting order parameter. Our finding may be best summarized in Fig. 13. 

FIG. 13. (a) The best scan path to test our effects is the path denoted by "scan l."(b) The expected LDOS profile measured 
on the bond-center points (O sites on a Cu02 plane) around a single SU(2) vortex with the real band structure of BSCCO 
being taken into account. For an illustration we assumed £c = S and Xa ~ 10. The staggered peak-dip (SPD) structure appears 

over Ic. ~ r ~ \a and vanishes deep inside the d-wave SC state (r 3> Xa)- We expect almost no effects along the scan path 
denoted by "scan 2." 

The best scan path to test our effects is shown in Fig. 13(a) as "scan 1, "along which the LDOS on the bonds exhibits 
specific peak and dip structure alternating from bond to bond in a staggered manner as indicated in Fig. 13(b) [the 
LDOS shown here is obtained under the same setting as in Fig. 9]. The SPD structure appears over the region 
^ ^ Aq and vanishes deep inside the rf-wave SC state. The core size tc presumably extends over a fermion 
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coherence length S^p vf/Aq which may amount to a few lattice scalesc^ and the energetics of a single vortex 
supports that large value of Xa/^c tends to be favored. Thus we are hopeful that there is certainly the region 
^ r ^ Xa over which our effects are detectable. Due to the lattice symmetry, the unit cell doubling effects on the 
LDOS is detectable only on the bonds. Thus, we have just a typical V-shaped profile of bulk d-wave SC along the 
path denoted by "scan 2" in Fig. 13(a). Although the qualitative feature of the LDOS profile may not be so sensitive 
to the doping x in the underdoped regime, fine detail of the SPD structure depends on the doping dependence x which 
controls gq oc x and Aq. In particlular, existence or absence of the second peak/dip depends on x. Nevertheless, 
we have at least one resonance (or a pair of van-Hove singularities) on the reduced zone boundary at the energy 
around lo = ±(2xo>/i3 + oo) [Eq- ( 4.30| )] which always lies inside the V-shaped profile for reasonable band structure 



parameters of BSCCO. 

As for an experimental setup for BSCCO sample, the best place to test our prediction is the O site around the 
vortex center on the Cu02 plane. The size of the Wannier function at the O sites on Cu-O-Cu bonds is presumably an 
order of 10% in magnitude as compared with the nearest Cu sites. However, it is noteworthy that the STM tunneling 
into the O sites may take place directly via the STM tips, while the tunneling into the Cu sites on the Cu02 plane 
takes place indirectly through the Bi atom on the BiO layer .0 Thus, we are hopeful that the STM signal may more 
sensitively detect the LDOS profile at the O sites than at the Cu sites. We stress that the SPD structure is totally 
ascribable to the unit cell doubling and the robust topological texture. Therefore, we may safely say that the SPD 
structure survives any tunneling matrix element effects and can directly be detected through STM experiment. 
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APPENDIX A: INTERNAL PHASE OF THE LOCAL BOSON CONDENSATE 

The uniform d-wave SC state in the d-wave gauge is described by, 

Uf^=-Xor' + {~iy-+'-Aon, (Al) 

hm=^'f), (A2) 

which is just equivalent to the U(l) MF solution for the d-wave SC state. Now, thanks to the SU(2) symmetry, the 



same state can be described in the SF gauge via the SU(2) gauge transformation w] given by Eq. ( 2.10 ). The gauge 
transformation converts Ufj and hoi to 

U^/ wjuf^wj = -At"^ exp [z(-l)'-+J«$or3] , (A3) 

where A = a/xo"+A^, and $o = tan^^(Ao/xo)- Now, the low energy excitations around the SC state in the SF gauge 
are obtained by fixing Uff and then rotating the boson condensate in the internal SU(2) space. The direction in the 



internal SU(2) space is specified by the internal angles and 9 as in Eq. ( 2.12 ). We obtain this parameterization 
more directly through transforming h^^ by 
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i^+iy 1 





3 


exp 





gl = exp 

which we will encounter in Sec. Ill A [see Eq. (|3.l| )] when we make a gauge transformation. 

APPENDIX B: LONDON MODEL ANALYSIS OF A SINGLE SU(2) VORTEX 



(A5) 



Here we apply the London model prescriptiorH to a single SU(2) vortex. Plugging Eq. ( ^.15 ) into Eq. ( 2.16| ) gives 
Fk ^ Fv + Fy, where 



Fv^ 
Fy = 



X 
X 

2mh 



drV{r)^, 



with 



V = -V(y59 A = Va A, 

2 c c 

1 o 1 , 

2 2 



(Bl) 
(B2) 

(B3) 
(B4) 



being the superfluid velocities associated with A and fields, respectively. This decomposition indicates that A and 
gauge fields are decoupled at the mean field level. The stationality condition with respect to A, 5{Fv + Fa)/5A — 0, 
is reduced to ^^H{r) — X^'^H[r) = — ■^j5-e2(5(r), which gives the solution 



H{r) = 



271 Xi \Xi 



(B5) 



where Kq is the zero-order modified Bessel function of an imaginary argument. The London penetration depth is 
defined as A| = rniiC^ /Ane^x. The EM unit flux is (/i™ — hc/2e where we retreived Plank constant. The physical 
supercurrent associated with A becomes 



r m r 



r 



(B6) 



which globally circulate around the vortex center over the length scale . 



Taking account of Eq. (2.21), the stationality condition with respect to a^, 6{Fy + Fa)/Sa^ — 0, is reduced to 
-2xtv{r) + (J J dr''VrK{r — r') x h{r') = 0. Taking curl of this equation and going to Fourier space, we obtain 



h{q) 



/.gauge 



(B7) 



'l + XaHqq^ l + Aalql' 

where we made use of Eq. ( ^.221 ) . The gauge field penetration depth and the unit flux associated with it are given by 



Aq = TTifjCr/a;, and (f)^^^'^ ~ ft,/2, respectively. Fouier transform of Eq. (B7) gives 

h{r) 



^gaugc 



XI 



-Nnl — 

Xn 



(B8) 



where Ho(z) and No(2;) denote the Struve functiorO and the Bessel function of the second kind, respectively. We 
note that Aa/A^ ~ ^ 1, where we used (t ~ \/ J A ^ l/^/rub (Al reaches ^ 500 with lattice unit in BSCCO). 
The essential point here is that we can reasonably assume the EM gauge field H extends much broader than the 
internal gauge field h. Under this circumstance, the effect of the EM gauge potential A is negligible as compared 
with the internal gauge potential a"^. That is to say, in our vortex model, a half flux quantum of the EM gauge 
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field, v4, penetrates over a huge region r < A^,, as compared with a half flux quantum of the internal gauge field, a^, 
which penetrates over a region r < Aa- In Fig. 14, we show spatial decay of H {r) / {(j)'^ /2t:\\) and h{r) / {(j)^^^'^ / \^^) 
assuming Al/Aq = 50 and Aa = 10 with lattice unit. The apparent divergence of h{r) at r = should be taken as an 
artifact of continuum limit, since there is natural cutoff of an order of the inverse lattice scale. 

FIG. 14. Spatial dependence of H{r) / {4)'^^ / 2-^X1) and h{r)/{<j>f'^'^° assuming Al/A„ = 50 and Aa = 10 with lattice 
unit. Apparent divergence of h(r) at r = should be taken as an artifact of continuum limit. 

We shall now argue how A^ can become larger than ic in terms of energetics of a single vortex. The energy associated 
with a single vortex consists of the following contribution: the cost for the SF core formation Ccore, the electromagnetic 
contribution cem = Fy + -Fa, and the contribution of the internal gauge field Cgaugc = -Fu + ^a- The enegy cost for the 
SF core formation is estimated as as0 ecoro ^ V J /S.H^x^l'^ , which favors smaller ic- On the other hand, the core size 
cannot be smaller than x^^/"^ without costing too much kinetic energy. Thus, we conclude that the SF core occupies 
a radius of Ic ^ x^^l"^ at MF level. In the present scheme, it is quite reasonable to expect that as the doping x 
decreases, the core size becomes larger because the energy difference between the SC and the SF state decreases as 
a; ^ 0. 

The electromagnetic contribution comes from eEM = Fv + Fa which reduces toS 



EEM = ^ / dr + Ai(V X Hf\ 



4m{, I, 



In^. (B9) 



To compute Cgaugo we first take curl of Eq. (B4) and go to Fourier space to obtain 



and similarly 



Recalling x/2ra}j ~ cr/2Aa, we have 



v{q) = -nT —-——^ BIO 

q- 1 + Xaq 



\ . (Bll) 
q^ 1 + Xaq 



gauge 



El 1 TTC An 
^ — = ^lni;^. B12 
ql + Xaq 2Xa 4 ^ ' 



We thus have the energy cost associated with a single vortex: 



^vortex — ^corc ~r f^gaugc £EM 



This result is consistent with a little bit more qualitative discussionjs^ i.e., a standard hc/2e vortex is possible with 
the SF core which does not cost too much energy as x ^ 0. The SF core size £c would like to be small as possible 
with the lower bound ic ^ x^^^^, while the size of the gauge field distribution Aa = nibcr/x would like to be large. 



APPENDIX C: COMPUTATION OF LDOS 



1. Derivation of Eq. (4.14) 



We start with Eq. ( |4.3| ). The LDOS at the plaquette-center, site-top, and bond-center points, as indicated in Fig. 6, 
are given by 
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2x 



i j"=l,2,4,5 



TT 



iV''™^(rB,u;) 



(CI) 

(C2) 
(C3) 



respectively. The envelope function is simulated by ai{r) = e I'' = 1/2 with lattice unit) in the bond 

direction. Tuning ai(ri) = 1 at the site-top points, we put Q!i(t"c) — a2{rc) — a4{rc) — cisifc) ^ = ac, and 



ai(^B) ~ oi2irB) = c^B- Using Eq. (4.11), we obtain 

4 



^Al + l/'2 + V'4 + V'5 
^1 



V2 



feeRZ 



A/2; , , , Az-y 

cos — COS yVfc ± sm — sm —ipk+Q 



fceRZ 



COS Y'V^fe ± * sin ^ipk+Q 



(C4) 
(C5) 
(C6) 



where — and + signs being for the cases where the site 1 belongs to the A and B sites, respectively. Constructing the 



propagators with these wave functions, we obtain a formula Eq. (4.14) in a concrete form. For example, in Eq. ( |C3| ) 
we obtain 



(C7) 



The uniform and staggered counterparts are computed as 



g^{rB,iuj) = E E ^ 



i{k~k')-rB 



fceRZ fc'eRZ 



cos ^ cos ^Q^{k, k' , iuj) + sin ^ sin ^Q^{k + Q,k' + Q, iuj) 



(C8) 



and 



fceRZ fc'GRZ 



COS y sin ^g^{k, k' + Q, lio) - sin y cos ^^^(fe + Q, k' , iu) 



= Y • t-b) sin 

fcGRZ q~small 



— ^ COS — ^ 



gF(fc_|,fc + | + Q,zc.) + gP(fc + | + Q,fc-|,zc^) 



(C9) 



which give Eq. (4.15) and Eq. (4.20), respectively. At the last step we replaced k — k' and (A; + k')/2 with q and fc, 
respectively. 



2. Computation of the generalized propagators 



Next we compute 



g'{k + ^ + Q.k- |,r) = -(T,Vfe+^+Q(r)V4_,), 
S'^ik - |, fc + I + Q, r) = -(TW^fc-i (r)^I+^+Q), 
at Born level. First, we rewrite the perturbation term Eq. ([l.7|) in momentum space: 



(CIO) 
(Cll) 
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= H ^i'^ 

feeRZ fe'GRZ 

+ iJ ^ ^ sin 
fceRZ fc'eRZ 



(C12) 



2 



y'^'^y \ vy{k - fe') Ui(AoT' - XoT^)^k'+Q - V4+q(Aot' - XOT^)^k' 



By replacing k ~ k' and (fc + A;')/2 with q and A;, respectively, and recalling Eq ( B1C| ), we reach Eq. (4.12). The first 
order contribution of 6H^ to the propagator is obtained as 



+ f + Q, fc - f , r) = - ^ e^-^g^ik + I + Q, ■iL,)Ck+Q{q)g!{k -^,^u;), 



(C13) 



I.e. 



g^{k + ^+Q,k-^,iuj) = -g^{k + I + Q, iu;)Ck{q)g^ik 



Similarly, we obtain 



- |, fc + I + Q, r) = -g^{k - ^,tu;)Ck{q)g^{k + | + Q, i^)- 



Recalling Eq. (4.8), we explicitly write down perturbative corrections: 



gF^k - |, fe + I + Q,zc.) + gF(fe + I + fc _ |,,^) 



1 



1 



1 1 





+ E+ 


iuj 


+ E- 




1 




1 






iuj 


+ 




1 




1 


ILU 


+ E+ 


iuj 





- [u+Ck{q)U- + f/-Cfc(<j');7-t 
[y+Cfe(q)T/_ + x/_Cfe(q)y+] 

[C/+Cfc(q-)V^_+y_Cfc(q)C/+] 
[V+Ck{q)U+ + U.Ck{q)V+], 



where 



E+ — El 
U± = \ + 1±T^ + ri±r^] 



k-i 



(C14) 



(C15) 



(C16) 



(C17) 



(C18) 



with 7+ = -ik+q/2+Q/Ek+q/2+Q, 1- = lk-q/2/ Ek-q/2, V+ = Vk+q/2+Q / Ek+q/2+Q, and T]- = Vk-q/2 / Ek-q/2- Taking 
11 component of Eq. ( pl6 ) and then performing analytic continuation, iuj ^ lu + iS, we reach Eq. ( |4.24D . 



APPENDIX D: EXPLICIT FORM OF G"" 



By simply taking inverse of the matrix iujl — with given by Eq. (4.34), we obtain an explicit form of 



G^{k,iuj). The 11 and 33 components are given by 



[G^{k,iu;)]n = [G(fc + Q,zc^)]33 
(ioj)^ + Ak{iuj)^ — BkioJ - Ck 



D{k,iLu) 
Uik , Vik 



U2k 



V2k 



iuJ - -Efc iuj + E^. iuj — E^ iuj + E'^ ' 



(Dl) 
(D2) 
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where D{k,iuj) = det[iu) — T^] = {iuj — E^, ){iLu + E^, ){iu) — E'^){iuj + E'^), and 
Ak ^ ak- Ik, 

Bk = K+7fe)'+^/I + A| + Ml, 

Ck = (flfe + 7fe)(a| - 7fe - A| + /Ufc) + (flfc - 7fc)?7^ + 2r]kXkfJ.k- 
The generahzed coherence factors are given by 

Uik = 2^ [i^k)'^ + ^k^k - Bik - Cik/E^] , 
^1'= = ^ [(^fc )' - Ej^^k -Bk + Ck/E^] , 
U2k = \jr^ HKf - EtAk + Bk + Ck/E+] , 
V^k^l^ [-(Eif + E+Ak + Bk- Ck/E+] , 

where 



Pfc ^ [E-f {E+f = -AJali^l + XI) + ivkXk + Ik^ikV. 



(D3) 
(D4) 
(D5) 



(D6) 
(D7) 
(D8) 
(D9) 

(DIO) 



Similarly, 



where 



[GF(fc,zc^)]i3 = [GP(fc,ia; 

{iu)'^ + 2akiuj + Ck 
D{k, luj) 

Uik , Vik 



J31 



= -iXk 



U2k 



V2k 



iuj - E. 



iu} + E 



k 



iuj - El 



iu} + E2 



(Dll) 



Ck=al--fl- Xl+ r/l - + 27fer/fe/ife/Afc, 



(D12) 



Vik^-- — 



1 1 
2Pk 

V2k--'-^ 

2Pk 



E^ + 2ak + Ck/E^ 
^Ef, + 2ak — CkjEy, 
-E+ - 2ak - Ck/E+ 
E+ -2ak + Ck/E+ 



Now, LDOS at the midpoint on the bond connecting and + is given in the form of Eq. (4.38) with 



7Vo(w) = ^cos2 ^ [Uik5{uo - E^) + Vik5{uJ + ) + U2k5{u^ - E+) + V2kS{u; + E+)] , 

k 

SNiij) = ^ Afc sinfc^ \UikS{u; - E^ ) + Vik5{Lo + E^ ) + U2k5{'^ - E+) + V2k5{Lo + E+) 



(D13) 
(D14) 
(D15) 
(D16) 

(D17) 
(D18) 



fceRZ 



These equations further reduce to Eqs. ( ff.44 ) and ( 4.45 ) 



26 



^ p. W. Anderson, Science 235, 1196 (1987). 

^ I. Affleck and J. B. Marston, Phys. Rev. B 37, 3774 (1988). 

3 G. Kotoliar and J. Liu, Phys. Rev. B 38, 5142 (1988). 

* Y. Suzumura, Y. Hasegawa, and H. Fukuyama, J. Phys. Soc. Jpn. 57, 2768 (1988). 
^ G. Baskaran and P. W. Anderson, Phys. Rev. B 37, 580 (1988). 

® L. B. loffe and A. I. Larkin, Phys. Rev. B 39, 8988 (1989). 
P. A. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1992). 

* L Affleck, Z. Zou, T. Hsu, and P. W. Anderson, Phys. Rev. B. 38, 745 (1988). 
" D. Poilblanc and Y. Hasegawa, Phys. Rev. B. 41, 6989 (1990). 

° M. U. Ubbens and P. A. Lee, Phys. Rev. B. 46, 8434 (1992). 
1 H. Fukuyama, Prog. Theor. Phys. Suppl. 108, 287 (1992). 
^ X.-G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). 

^ P. A. Lee, N. Nagaosa, T. K. Ng, and X.-G. Wen, Phys. Rev. B 57, 6003 (1998). 
■ D. A. Ivanov, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. 84, 3958 (2000). 
P. W. Leung, Phys. Rev. B. 62, R6112 (1988). 
^'^ S. H. Pan, E. W. Hudson, and J. C. Davis, Rev. Sci. Instrum. 70, 1459 (1999). 
^''L Maggio-Aprile et al., Phys. Rev. Lett. 75, 2754 (1995). 

Ch. Renner, B. Revaz, K. Kadowaki, L Maggio-Aprile, and O. Fischer, Phys. Rev. Lett. 80, 3606 (1998). 
1^ S. H. Pan, E. W. Hudson, A. K. Gupta, K.-W. Ng, H. Eisaki, S. Uchida, and J. C. Davis, Phys. Rev. Lett. 85, 1536 (2000). 
" N. Nagaosa and P. A. Lee, Phys. Rev. B 45, 966 (1992). 
S. Sachdev, Phys. Rev. B 45, 389 (1992). 

A. Himeda, M. Ogata, Y. Tanaka, and S. Kashiwaya, J. Phys. Soc. Jpn. 66, 3367 (1997). 
H. Tsuchiura, Y. Tanaka, M. Ogata, and S. Kashiwaya, J. Phys. Soc. Jpn. 68, 2510 (1999). 
M. Franz and Z. Tesanovic, Phys. Rev. B 63, 064516 (200 1). 
2^ J. H. Han, Q.-H. Wang, and D.-H. Lee, | cond-mat001245C . 
Q.-H. Wang, J. H. Han, ar id D.-H. Lee, cond-matO 10204^ . 
P. A. Lee and X.-G. Wen, 2ond-ma t/0008419| . 

T. Hsu, J. B. Marston, and L Affleck, Phys. Rev. B. 43, 2866 (1991). 
J. Kishine, P. A. Lee and X.-G. Wen, Phys. Rev. Lett. 86, 5365 (2001). 
3° E. Dagotto, E. Fradkin, and A. Moreo, Phys. Rev. B 38, 2926 (1998). 
X.-G. Wen, Phys. Rev. B 44, 2664 (1991). 



15 



31 



' When we write Eq. (2.9) as Uff = [ gp ] , an additional gauge transformation xW ~^ *xff ^^'^ xW* ~^ ~*xff * 

y Xtj J " 

leads to an alternative gauge choice: Uff — —ixo— (— l)'°'''"-'''r^Ao, which was used in the papers by Lee-Nagaosa-Ng-Wenjlj 
and Lee- Wen £J This gauge transformation just causes (7r/2, 7r/2) shift of the Borillouin zone and leads to the 6i and 62 boson 

bands having the same minimum at (71/2, 7r/2). Then, the local boson condensate is written as hf^ = y^e"*'*^^'"''^'^^ ( 



instead of Eq. ( I2.I4 ) . This change of the gauge choice is of course just apparent and does not change any physical consequence 
of our vortex model. 
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The following formula is useful to evaluate the Fourier transform of Eq. (|B7|): xJo(ax) _ i _ l[Ho(a) — No(a)] for 

a > 0, where Ho{z) and No(z) denote the Struve function and the Bessel function of the second kind, respectively. 

To have some feeling of this situation, we consider a non-interacting electrons on a linear chain described by a toy Hamiltonian 

H = ^X/J* ^~ (^l)*<5] [cj+ic''"'^-' '"^Ci 4- H.c.]. It is easy to see that the staggered orbital currents flow on bonds, jij = 

it{c'.^Cia — c\^Cjcr), if 5 = aud 7^ 0, while the LDOS on the bonds may become staggering if 5 7^ and 4> — Q. The orbital 

current flowing and the staggered L DOS are complementary to each other. 
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